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Abstract 

In this work we analyze a system consisting in two-dimensional position-dependent massive 
particles in the presence of a Morse-like potential in two spatial dimensions. We obtain the exact 
wavefunctions and energies for a complete set of eigenstates for a given dependence of the mass 
with the spatial variables. Furthermore, we argue that this scenario can be play an important role 
to construct more realistic ones by using their solution in perturbative approaches. 
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I. INTRODUCTION 


Some years ago, the materials science took an important step forward by obtaining the 
fabrication of small conducting devices known as quantum dots (QDs) [1]. In those devices, 
it is possible to conhne several thousand electrons in a small region whose linear size is about 
0.1 — 1pm [2], The fundamental characteristic of QDs is that they are typically formed by a 
two-dimensional electron gas, where, by applying an electrostatic potential, the electrons are 
conhned to a small region, which is called “dot”, in the interface region of a semiconductor. 
A very important advantage of QDs is that their transport properties are readily measured, 
allowing an experimental control. Moreover, the effects of time-reversal symmetry breaking 
can be easily measured by applying a magnetic held [3]. Nowadays, a variety of theoretical 
and experimental research about small conducting devices, such as QDs, has been the focus 
of many scientists and engineers attention [4-7]. From a phenomenological viewpoint, QDs 
are very small structures, where the laws of quantum mechanics (QM) are the most impor¬ 
tant ingredients to describe their properties. Thus, as a natural consequence of practical 
applicability of the theoretical framework of QM, a great interest arises for exact solutions 
of two-dimensional conhned systems, which can be fundamental to explore the physics in 
small conducting devices, such as the QDs. In the light of these facts, it was shown in Ref. 
[8] that it is possible to hnd exact solutions of the two-dimensional Schrodinger equation 
with the position-dependent mass (PDM) for the square well potential in the semiconductor 
quantum dots (SQDs) system. Another important work in this context, it was presented 
by Schmidt, Azeredo, and Gusso [9], where the authors have studied both the problems of 
quantum wave packet revivals on two-dimensional inhnite circular quantum wells (CQWs) 
and circular quantum dots (CQDs) with PDM, showing the results for the eigenfunctions, 
eigenenergies and the revival time for spatially localized electronic Gaussian wave packets. 
At this point, it is important to highlight that the importance in adding a PDM is due to 
the fact that the system will take into account the spatial variation of the semiconductor 
[10-17]. However, as consequence of inclusion a PDM, the system becomes ambiguous at 
the quantum level, and the ordering ambiguity problem (GAP) is one of the long standing 
unsolved questions in quantum mechanics. As we know the GAP has attracted the attention 
of some of the founders of the quantum mechanics, namely. Born, Jordan, Weyl, Dirac and 
von Newmann worked on this problem, as can be verihed from the review by Shewed [18]. 
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This is viewed as a deep problem in QM, which has advanced very few along the last decades. 
But all is not lost, it was shown that the ordering ambiguous problem has a very special 
importance for the modeling of some experimental situations like electrons in perturbed 
periodic lattices [19], impurities states and cyclotron resonance in semiconductors [20], the 
structure of electronic excitation levels in insulating crystals [21], the dependence of nuclear 
forces on the relative velocity of the two nucleons [22, 23], and more recently the study of 
semiconductor heterostructures [12, 25-28]. Moreover, some time ago, it was discussed in 
the literature the exact solvability of some classes of one-dimensional Hamiltonians, where 
the potentials has a PDM, with ordering ambiguity [29], after that, a large number of works 
regarding one-dimensional Hamiltonians with ordering ambiguity has emerged in the scien¬ 
tific community along the last few years [30-36] . Another interesting research line regards to 
the supersymmetry approach to one-dimensional quantum systems with spatially-dependent 
mass, by including their ordering ambiguities dependence [37-61]. On the other hand, as 
far as we know, some physical systems like ones where a magnetic field is present [62-65], 
lead naturally to the necessity of a two-dimensional analysis. In the face of this situation, 
it was presented in Ref. [66] a general approach for the problem of a particle with PDM 
interacting with a two-dimensional potential well with finite depth, where the ordering am¬ 
biguity was taken in account. In that work, it was shown that the considered system retain 
an inhnite set of quantum states, which usually do not happens in the case of the constant 
mass systems. Furthermore, it was verified also that the SU(2) coherent state corresponds 
to a stationary state. Also recently, numerous other theoretical studies have been con¬ 
ducted on two-dimensional position-dependent mass Schrodinger equation (PDMSE). These 
include the two-dimensional quantum rotor with two effective masses [67], kinetic operator 
in cylindrical coordinates [68], exact solutions for the PDMSE in an annular billiard with 
impenetrable walls [69], and a particle with spin 1/2 moving in a plane [70]. 

Here, we will address the position-dependent mass (PDM) type of the systems in two 
spatial dimensions (2D) by using Cartesian coordinates. We will introduce a very interesting 
system where, as we are going to see below in the manuscript, the relation between the 
quantum numbers introduced along the procedure of resolving the equations of the system 
and the energy eigenstates organization is somewhat remarkable. 

This paper is organized as follows. In Section H, we review the effective Schrodinger 
equation in two-dimensional Cartesian coordinates. In Section HI, we introduce the Position- 
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dependent massive particle with Morse-like terms and its exact solntions. In Section IV, we 
present our conclusions and directions for future work. 


II. EFFECTIVE SCHROEDINGER EQUATION IN TWO-DIMENSIONAL 
CARTESIAN COORDINATES: A BRIEF REVIEW 


In this section, we will recapitulate the results presented some years ago in Ref. [17]. Let 
us start with the ordering dehned by von Roos [13, 29] for the Hamiltonian operator, which 
in one-dimensional space is written in the following form 

^ = 1 [M^pM^pM^ + + V{x), (1) 

where p is the momentum operator and M = M{x) is the position-dependent effective mass. 
Moreover, a, (3 and 7 are arbitrary ordering parameters which must to obey the relation 

a + -f 7 = - 1 . ( 2 ) 

At this point, it is important to highlight that the above relation is necessary to get the 
correct classical limit. 

Applying the canonical commutation relations, we have 


M' (M'V M” 

= — - ih{P + ‘^oc)—p - tfa{P + a- 1 )^^ - 

Through the relation (2), the effective Hamiltonian operator [29] is given by 


( 3 ) 


H = 


;p + 


ihM' 


2M^“ ' 2 M2 

where the effective potential U{a,'y,x) is written as 


p U{a,'y,x) -F V{x), 


( 4 ) 


f/(a, 7 ,x) = 


h2 


4M3 


(a -|- 7 )M 


d'^M 

dx'^ 


. , fdM 

2 (« + 7 + «7) ( ^ 


Therefore, we can now write the effective Schroedinger equation in the form 


( 5 ) 


h2 

2 M{x) dx^ 2 


dM/dx 

M2 


^ + [V{x) + U(a, 7 , x) - E]^ = 0. 


( 6 ) 


In the case of a set of two-dimensional Cartesian coordinates, where M = M{x,y), the 
effective Hamiltonian operator 
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H = 




• fe / (M.ri 4- dM^ 

* ^ / dx ^ dy Py 


M2 


+ U{a.-f,x) + V{x,y), (7) 


where U{a,'y,x,y) is the effective potential. Now it can be written in the form [17] 


U{a,-^,x,y) = <j (a + 7) 


MxX + Myy 

M 


2{a + 7 + 07 ) 


M ) ^ \ M 


with Mx = dM/dx and My = dM/dy. Therefore, we have 


, ( 8 ) 


^ ^ ^ x,y) + V(x, y), 


where, in this case 


U{a,-f,x,y) = 




4M 


V^M 

-2(a + 7 + a7) 




M 


M 


In the next step we can use a typical Schroedinger equation 

^‘^X + Veff{x,y)x = Ex, 


( 9 ) 


( 10 ) 


( 11 ) 


2M{x,y) 

and if x{^,y) = y) is the solution of it, the equation above can be rewritten as 

follows 


h2 






2M{x,y) M{x,y) 

V{x,y) 




VV + 


7 .)!} 


2M{x,y) 

The above equation have a Hamiltonian operator dehned by 


^|J = Eifj. 


( 12 ) 


H = 






2M{x,y) M{x,y) h 
Note that we can choose 

i 


(^(T).y + V 




2M{x,y) L 


VV + ('^a)2 




(13) 


( 14 ) 


such that 


5 


























(15) 


M 


= -2^^ 


a. 


Thus, we have 


a = ln{M 2). 

Now, we may rewrite the equation (13) as 
rr 1-^2 ih^ M . 

2M ^ 2M M ^1 AM 
On the other hand, the wavefunction is re-scaled as 
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2 \ ~W 


V^M 

M 


(16) 


(17) 


^jJ = Ma y. 


From these results, we see that [17] 


Veff{x,y) = V{x,y) + 


AM 


2(q!-F7 + Q!7-F- 




M 


M 


- (a -f 7 -h 1) 


V^M 

M 


(18) 


(19) 


Finally, we must comment that for an equivalent system with constant mass, the equation 
(11) can be written as 


TrV‘^X + U^,fX = ix, 


( 20 ) 


with ^ constant and 


Ueff-i = M{x,y)V{x,y) + — 


(a 7 -h 1) 


V^M' 

M 


2\a + ^ + a^ + ^ 


EM{x,y), 




M 


M 


+ 


( 21 ) 


Through the above result, it was studied in [17] the problem of a particle with a position- 
dependent mass interacting with two-dimensional potential well with hnite depth, as well as 
under the influence of a uniform magnetic field. There, it was discovered that the system 
retains an inhnite set of quantum states. In the next section, we explore the problem where 
the PDM is Morse-like. 
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III. POSITION-DEPENDENT MASSIVE PARTICLE WITH MORSE-LIKE 
TERMS 

An important problem in quantum mechanics is that one related to the vibrations of 
diatomic molecules, and the case of vibrations of a two-atomic molecule are well described 
by the Morse potential. On the other hand, there is a growing number of applications of 
quantum wells and quantum dots. In fact, those systems present a small spatial region 
capable to conhne quantum particles. As one can see in Figure 1, the bidimensional Morse- 
like potential can simulate such kind of physical situation and it has the advantage, as we 
will see below, of being exactly solved. Furthermore, having the exact solutions in hands, one 
can use them in order to describe more realistic problems by using approximation techniques 
which make use of those exact solutions. Therefore, with this motivation in our mind, in 
this section, let us present an example which can be exactly solved. Thus, we will consider 
that 


M(x, y) = M„ [l + gie-"'* + gae-”® + gze-"”" + 


( 22 ) 


Note that the spatial dependence of the mass is similar to that of a Morse potential in 
two dimensions. Thus, plugging this mass in the formula given by (21), we obtain 

Ueff-^ = Mo[l+ + g^e-^^y + + g,e-^^^y] V{x, y) 

-EMo [1 + + g^e-^^y + 


( f 3 

2 I a + 7 + a7 + - 


[algfe -1- a^g^e _|_ 4a‘lg2e -|- 4Q;2fi'4e 

(1 + + gse-^^^y + + g^e-'^^^yy 


-( q ; + 7 + 1 ) 


afgie -|- a^g^e + 4:alg2e + 4Q;2fi'4e 


1 + gie-°‘^^ + g-ie-°‘^y + g2e-‘^°^^^ + g^e-'^^^y 
In order to work with a exactly solvable model we can assume the following ordering 


(23) 


a + 7 + 1-0, a + 7 + Q:7+--0, 


(24) 


7 



whose solution is given by 


1 1 

a = —, 7 = —, 3 = 0. 
2 2 


( 25 ) 


In this way, we then obtain 


^2 


1 1 .y^2 1 


2M Tm' 

Consequently, in this ordering, the effective potential is written as 


(26) 


= M„ [1 + + 936 -““ + 926-"““ + 946 -""“] V(x, y) - EM„ 


£:A/„(9i6-““ + 936-"“ + 926-"““ + 946-""“). 


( 27 ) 


As an example, we can choose a potential under which the particle with position- 
dependent mass is moving. Then, here we will work with the following potential 

A + -f 53e-“22/ + 526-2“!^ -f 


V{x,y) = R + 


(28) 


Mo (1 -h -f gse-^^y + 5(26-2“!^ -f g 4 e-‘^°^^y) ’ 

where i? is a constant. In Figure 1, we plot a typical case where this potential can conhne 
particles. Therefore, through this potential, we can obtain the effective potential 


U,ff-^ = A + MoiR -E) + [B, + MoiR - E)g4] + [B^ + M^iR - E)g 2 ] e 


—2cx.\x 


+ [2B3 + Mo(i? - E)g^] e-^^y + [B4 + M^{R - E)g4] 
where we can easily see that ^ = —A + Mq{E — R). So that 


( 29 ) 


Ueff = Tl6-"‘- + 726-"““ + 736-"“ + 746-"““, 


with 


'yi = Bi + Mo{R - E)gi, i = 1,2, 3,4., 

Now, the Schroedinger equation (20) takes the form 

- V^x + ^ .y^e-^<^2y^x = £ X- 


( 30 ) 

( 31 ) 


( 32 ) 



where e = 

In order to solve the above equation, we can use the usual procedure of variable separation 


x{x,y) = X{x)Y{y). 
Thus, we get the equations for X(x) and Y(y) below 


(33) 


where 


dx^ 

^mX (x), 

(34) 

+ (h2e-“^" + U 2 e-^^^y)Y{y) 
dy^ 

= enY{y). 

(35) 

271 272 273 

_ 274 

(36) 

01 = . "1 = . 'fe = >'2 

“ ■ 


Furthermore, the energy spectrum is given by 

^mn T ^n- 


(37) 


Let us now determine the solution of X{x). Note that the equation (35) have the same 
form of (34), of course, written in terms of the variable y. In this way, it is necessary to 
solve only (34). Thus, we dehne the variable 2 : and constants y and A as 


c- ? • 




, A: = 


hi 


(38) 


cq ' tti 2ai\J\vi\ 

with -00 < X < 00. In this case, bound states are possible only for z/i > 0 and 171 < 0. Then, 
we have 

1 2 

e™ =-—[|? 7 i| - aivdq( 2 m + 1 )] , with m = 0,1, 2, 3,..., rrimax- (39) 


Furthermore, the function X{x) is given by 


X-miX) 




exp 




V\^\ - 


(y.\x 


(y.\ 


z -» 


2v^ _ 


aix 


ai 


( 40 ) 


where C‘^{x) are the Laguerre polynomials. Here, it is important to remark that the number 
of discrete levels is hnite and determined by the condition 


Ihil > al^/iy{2m max + !)• 


(41) 
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This happens due to the fact that the potential goes asymptotically to zero when x —)■ oo 
and its minimum value is negative. So the energy levels for bounded particles must be lower 
than zero, which leads to the above constraint. On the other hand, dehning 

/J • C- ^ ^ /\ • , ^ 

tti o>2 2q;2a/| zZ2| 

and solving the equation (35), we obtain 

[|? 72 | - a2v^(2^^ + 1)]^ , with n = 0,1, 2, 3,..., rimx- (42) 

4Z/2 

and 


Ym{y) 



with the condition 


fia2y + 


V\^\ 


a2y 


Oi2 




z -)■ 


2vT^ - 


022/ 


0.2 


\rj2\ > q;2v^( 2^ max + !)• 


Therefore, the total energy is written as 

£mn = {i£2 [Ihi| - aiy^(2m + 1)]^ + z /1 [|? 72 | - a 2 v/^( 2 n + 1)]^| . 

4Z/]^Z/2 ^ ^ 

Moreover, we write 


(43) 


(44) 


(45) 


Xmn{Xj y) 



(2VW\Y 

\ C.2 ) 


jjiOix T 



fj,02y + 


a/m , 


a2y 


0^2 




2a/N ^-oix^ £2/1 f 2a/H ^ 


—ct2y 


(46) 


Ol J \ 

We know that Emn = and ^mn = —A + Mo{Emn — R)- Consequently the energy 

eigenvalues will rise as solutions of the following transcendental equation 


8 ^2 i^mn) b4 )iA+ 

^mn) — b4 (^mn) 


lbs {^mn) \ — Oi\/ ^ ^ -^(2n + 1) 


+72 i^mn) 


lbs {^mn) \ — 02\j ^ ^ -^(2m + 1) 


(47) 
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where we defined Gmn = Mo{R — Enm), ai = h ai and a 2 = h 0 : 2 - Considering a 
symmetrical (in x and y) case of the mass and potential dependencies, in order to have a 
concrete example to study, we choose the parameters as given by: i? = 0, A = Mq = gi = 
93 = 0:1 = a 2 = H = 1, g 2 = gi = 0, Bi = B 3 = —1, B 2 = B^ = 1/8. In this case, the 
allowed energy levels are given in the Table below (note that due to the symmetry of the 
system, the pair (n, m) have the same energy as the one (m, n)). Furthermore, the potential 
profile appears in the Figure 1 and a plot where the energy spectrum is presented in scale 
appears in the Figure 2. Note that, since the energy of the bound state can not be lower 
than que smallest value of the potential and that this potential becomes asymptotically 
constant, for the case of the above parameters, the allowed values of the bound states shall 
be in the interval —0.40693 < En,m < 1- By observing both the Table I and the Figure 2, 


n 

m 

Ejim 

n 

m 

Efim 

n 

m 

Enm 

0 

0 

-0.0669873 

1 

3 

-0.161438 

3 

4 

0.250000 

0 

1 

0.250000 

1 

4 

0.250000 

3 

5 

0.531754 

0 

2 

0.433013 

2 

2 

-0.329156 

3 

6 

0.883975 

0 

3 

0.250000 

2 

3 

-0.116025 

4 

4 

0.410275 

0 

4 

-0.0188424 

2 

4 

0.170844 

4 

5 

0.631966 

1 

1 

0.661438 

2 

5 

0.542893 

4 

6 

0.910275 

1 

2 

0.957107 

3 

3 

0.0317542 

5 

5 

0.801042 


Table I: Energy levels. 

one can note that there are some interesting results in the spectrum. First of all, we observe 
that this potential presents a finite number of allowed bound states, which is not a surprise, 
since this already happens in the case of the one-dimensional Morse potential (even in the 
case with position-dependent masses). However, in the case analyzed, there are inversions 
of energies where sates labeled with higher quantum numbers present lower energies than 
states with lower quantum numbers, as happens in the case of atoms with somewhat great 
atomic numbers. On the other hand, beyond the some expected degeneracies, we observe 
that there is a eight-fold degenerated stated (the seventh exited one). In this case we checked 
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FIG. 1: Morse-like potential in two dimensions. 



FIG. 2: Energy spectrum. 


that one shall have an accidental degeneracy, since we checked that changing slightly some 
potential parameters this degeneracy disappears, becoming a quasi-degeneracy. 
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IV. FINAL COMMENTS 


In this work, we present a general construction of a class of a two-dimensional PDM 
systems in Cartesian coordinates, analyzing an exactly solvable case and discussing its or¬ 
dering ambiguity and some of their properties. We extend the idea to the problem where 
ones deal with increasing mass Morse-like. In this case we obtain the exact wave-functions 
and energies for a complete set of eigenstates. Since the energy of the bound states come 
from a transcendental equation, involving the quantum numbers of a pair of one-dimensional 
equations, we discovered that this system presents a behavior which emulates the inversion 
of excited states usually seen in atoms with high atomic numbers. Moreover, an interesting 
accidental degeneracy appeared. Finally, it is important to remark that one could use this 
exactly solvable system in order to construct more realistic ones by using their solution in 
perturbative approaches. 


Acknowledgments 

RACC thanks to UNESP-Campus de Guaratingueta and CAPES for hnancial support. 
ASD thanks to CNPQ for partial hnancial support, and J.A.O. thanks to DFQ of UNESP, 
Campus de Guaratingueta, where this work was carried out. 


[1] Y. Alhassid, Rev. Mod. Phys. 72 (2000) 895. 

[2] M. A. Kastner, Rev. Mod. Phys. 64 (1992) 849. 

[3] A. R. Wright, M. Veldhorst, Phys. Rev. Lett. Ill (2013) 9, 096801. 

[4] G. Anatoly, Phys. Rev. B 91 (2015) 20, 205105. 

[5] Y. Li, A. Kundu, F. Zhong, and B. Seradjeh, Phys. Rev. 90 (2014) 12, 121401. 

[6] B. Wahlstrand, I. I. Yakimenko, K.-F. Berggren, Phys. Rev. E 89 (2014) 6, 062910. 

[7] P. Tighineanu, M. L. Andersen, A. S. Sorensen, S. Stobbe, and P. Lodahl, Phys. Rev. Lett. 
113 (2014) 043601. 

[8] M. Lozada-Cassou, S. H. Dong, and J. Yu, Phys. Lett. A 331 (2004) 45. 

[9] A. G. M. Schmidt, A. D. Azeredo, and A. Gusso, Phys. Lett. A 372 (2008) 2774. 


13 



[10] D. J. BenDaniel and C. B. Duke, Phys. Rev. B 152 (1966) 683. 

[11] T. Gora and F. Williams, Phys. Rev. 177 (1969) 1179. 

[12] G. Bastard. Phys. Rev. B 24 (1981) 5693. 

[13] O. Von Roos, Phys. Rev. B 27 (1983) 7547. 

[14] Q. G. Zhu and H. Kroemer, Phys. Rev. B 27 (1983) 3519. 

[15] T. L. Li and K. J. Kuhn, Phys. Rev. B 47 (1993) 12760. 

[16] F. S. A. Cavalcante, R. N. Costa Filho, J. Ribeiro Filho, C. A. S. de Almeida, and V. N. 
Freire, Phys. Rev. B 55 (1997) 1326. 

[17] A. de Souza Dutra and J. A. de Oliveira, J. Phys. A: Math Theor. 42 (2009) 025304. 

[18] J. R. Shewell, Am. J. Phys. 27 (1959) 16. 

[19] J. C. Slater, Phys. Rev. 76 (1949) 1592. 

[20] J. M. Luttinger and W. Kohn, Phys. Rev. 97 (1955) 869. 

[21] G. H. Wannier, Phys. Rev. 52 (1957) 191. 

[22] 6. Rojo and J. S. Levinger, Phys. Rev. 123 (1961) 2177. 

[23] M. Razavy, G. Field, and J. S. Levinger, Phys. Rev. 125 (1962) 269. 

[24] G. Bastard, Wave Mechanics Applied to Semiconductor Heterostructres, Les Editions de 
Physique, Les Ullis, 1992. 

[25] C. Weisbuch and B. Vinter, Quantum Semiconductor Heterostructures, Academic Press, New 
York, 1993. 

[26] C. C. Wu, J. Sun, F. J. Huang, Y. D. Li, W. M. Liu, Europhys. Lett. 104 (2013) 27004. 

[27] M. Benito, A. Gomez-Leon, V .M. Bastidas, T. Brandes, and G. Platero, Phys. Rev. B 90 
(2014) 20, 205127. 

[28] J. Gong and Qing-hai Wang, Phys. Rev. A 91 (2015) 4, 042135. 

[29] A. de Souza Dutra and C. A. S. de Almeida, Phys. Lett. A 275 (2000) 25. 

[30] C. Quesne and V. M. Tkachuk, J. Phys. A 37 (2004) 426. 

[31] J. F. Carinena, M. F. Ranada, and M. Santander, Ann. Phys. 322 (2007) 434. 

[32] A. Ganguly and L. M. Nieto, J. Phys. A 40 (2007) 7265. 

[33] S. Choi, K. M. Galdamez, and B. Sundaram, Phys. Lett. A 374 (2010) 3280. 

[34] A. Arda, R. Sever, and C. Tezcan, Phys. Scripta 79 (2009) 015006. 

[35] A. Ganguly and A. Das, J. Math. Phys. 55 (2014) 11, 112102. 

[36] A. de Souza Dutra, J. Phys. A 39 (2006) 203. 


14 



[37] A. de Souza Dutra, M. B. Hott, and C. A. S. Almeida, Europhys. Lett. 62 (2003) 8. 

[38] A. G. M. Schmidt, Phys. Lett. A 353 (2006) 459. 

[39] C. Quesne, J. Math. Phys. 49 (2008) 022106; J. Phys. A 40 (2007) 13107; Ann. Phys. 321 
(2006) 1221. 

[40] S. M. Ikhdair and R. Sever, J. Mol. Str. Theochem 885 (2008) 13. 

[41] A. Schulze-Halberg, Int. J. Mod. Phys. A 22 (2007) 1735; 21 (2006) 4853; 21 (2006) 1359. 

[42] J. F. Carinena, M. F. Ranada and M. Santander, Ann. Phys. 322 (2007) 2249. 

[43] A. Ganguly and L. M. Nieto, J. Phys. A 40 (2007) 7265. 

[44] B. Roy, Mod. Phys. Lett. B 20 (2006) 1033. 

[45] O. Mustafa and S. H. Mazharimousavi, J. Phys. A 41 (2008) 244020; 39 (2006) 10537. 

[46] T. Tanaka, J. Phys. A 39 (2006) 219. 

[47] R. Koc, M. Koca, and G. Shaninoglu, Fur. Phys. J. B 48 (2005) 583. 

[48] U. Gunther, F. Stefani, and M. Znojil, J. Math. Phys. 46 (2005) 063504. 

[49] G. Chen, Chin. Phys. 14 (2005) 460. 

[50] G. Chen and Z. D. Chen, Phys. Lett. A 331 (2004) 312. 

[51] A. A. Stahlhofen, J. Phys. A 37 (2004) 10129-10138, 

[52] B. Bagchi, P. Gorain, C. Quesne, and R. Roychoudhury, Mod. Phys. Lett. A 19 (2004) 2765- 
2775, 

[53] K. Bencheikh, K. Berkane, and S. Bouizane, J. Phys. A 37 (2004) 10719. 

[54] J. A. Yu and S. H. Dong, Phys. Lett. A 325 (2004) 194. 

[55] C. Quesne and V. M. Tkachuk, J. Phys. A 37 (2004) 4267. 

[56] Y. C. Ou, Z. Q. Gao, and Q. H. Shen, J. Phys. A 37 (2004) 4283. 

[57] R. Koc and H. Tutunculer, Annalen der Physik 12 (2003) 684. 

[58] A. D. Alhaidari, Phys. Rev. A 66 (2002) 042116. 

[59] B. Roy and P. Roy, J. Phys. A 35 (2002) 3961. 

[60] S. Ramgoolam, B. Spence, and S. Thomas, Nucl. Phys. B 703 (2005) 236. 

[61] A. de Souza Dutra, M. B. Hott, and V. G. C. S. dos Santos, Europhys. Lett. 71 (2005) 166. 

[62] B. K. Cheng and A. de Souza Dutra, Phys. Lett. A 123 (1987) 105. 

[63] A. de Souza Dutra and B. K. Cheng, Phys. Rev. A 39 (1989) 5897. 

[64] A. de Souza Dutra, C. F. de Souza, and L. C. de Albuquerque, Phys. Lett. A 156 (1991) 371. 

[65] M. S. Abdalla and J. R. Choi, Ann. Phys. 322 (2007) 2795. 


15 



[66] A. de Souza Dutra and J. A. de Oliveira, J. Phys. A: Math. Theor. 42 (2009) 025304. 

[67] A. G. M. Schmidt, J. Phys. A: Math. Theor. 42 (2009) 245304. 

[68] O. Mustafa, J. Phys. A: Math. Theor. 43 (2010) 385310. 

[69] A. G. M. Schmidt, Phys. A 391(2012) 3792. 

[70] A. G. M. Schmidt, L. Portugal, and A. L. de Jesus, J. Math. Phys. 56 (2015) 012107. 


16 



